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SUMMARY 


A stability and control analysis is made of a class of inertia-wheel - 
controlled earth satellites whose orientation is controlled relative to a set of 
orthogonal axes, two of which are the local earth vertical and the normal to the 
orbital plane. The NASA Nimbus meteorological satellite is a typical example. 

It is demonstrated for a restricted but nevertheless important sub -class of 
these satellites, defined by a particular inertia distribution, that the sta- 
bility about the local earth vertical (yaw) and the stability of the vehicle 
angular momentum are independent of each other. This fact is used to extract 
many fundamental dynamical properties of the system, in particular, the influ- 
ence, on stability, of the orientation and time constant of the yaw sensing gyro. 
Furthermore, the analysis gives some insight into the role of gravity torques 
and into the inherent difficulties involved in the removal of angular momentum 
from the vehicle. 

Performance data for a satellite similar to Nimbus, obtained from digital 
and analog studies, are used to supplement the stability analysis and to 
demonstrate transient responses following control commands. 


INTRODUCTION 


An important class of satellites is one whose mission requirements can be 
adequately met only by a configuration which is stabilized so that a line fixed 
in it lies along the local earth vertical or close to it. Such a satellite is 
usually termed an "earth -oriented satellite" or simply an "earth pointer." Per- 
haps the most widely discussed type of earth -oriented satellite control system 
is one which uses a combination of inertia wheels for momentum storage and mass 
expulsion for momentum removal. Research into the behavior of this type of 
system has been reported by Roberson in references 1 and 2 and Cannon in ref- 
erence 3 . The work of these investigators has yielded considerable insight into 
the general dynamical properties of earth-oriented satellites controlled in this 
way. 


The studies reported here are also concerned with the problems and behavior 
of earth -oriented satellites controlled by inertia wheels. The aim, however, is 
somewhat different from that of the previously mentioned authors in that a 
specific type of inertia-wheel control system is investigated in considerable 
detail. The primary aim of this control system is to provide full three -axis 
stabilization of the satellite relative to a set of orbit -fixed axes and in 
addition some measure of angular control (yaw control) about the earth -pointing 
axis. The essential features of the control system are that it has three inertia 
wheels with mutually orthogonal spin axes, each inertia wheel receiving signals 
from a separate attitude sensor. The sensors for the two axes located in the 
local horizontal plane (pitch and roll axes) are horizon scanners and for the 
earth -pointing axis (yaw axis) the sensor is a single -degree -of -freedom gyro. 

The report is divided into two major sections. The first deals with general 
theoretical considerations and the second with computed results, both digital and 
analog, obtained on a typical example of this type of control system. The numer- 
ical data used for the example were representative of the Nimbus meteorological 
satellite at one stage in its development. 

The theoretical section commences with considerations basic to a complete 
understanding of the behavior of the system. Included are some general results 
for very small and very large control system gains, a discussion of the best 
gyro orientation from a steady -state point of view and some steady -state results 
for large angles of yaw. They are followed by an analysis of the linearized 
equations describing the satellite system both from the stability and steady- 
state aspects. It is here that the main contributions of the report are to be 
found. The theoretical part of the report is completed with an analysis which 
demonstrates the feasibility of minimizing the effects of certain control system 
failures, in particular, a complete failure of the yaw channel. 

The section of the report containing the computed results on a typical 
satellite is organized so that the topics are introduced in roughly the same 
order as in the theoretical section. These results are used to strengthen the 
theoretical arguments and to reveal problems not readily apparent from the 
theoretical analysis . 


SYMBOLS 

Reference Frames 


All reference frames are right-hand systems. 

Oj_ orbital reference frame; 6 1; .along the velocity vector of the center 

of mass, o 2 normal to the orbital plane, o 3 along the local 
vertical, positive towards the center of the earth (See fig. 1.) 

bj_ vehicle principal axes of inertia; along roll, b^ along pitch, 

b 3 along yaw 


2 


m,* 


=xi 


a yi 


"marked" reference frame, oriented with respect to the b^_ frame 
through a constant small error angle matrix 

actual satellite fixed reference frame of the gyro, oriented with 
respect to the mq frame through a small error angle matrix 

actual satellite fixed reference frame of roll sensor, oriented 
with respect to the mq frame through a small error angle matrix 

actual satellite fixed reference frame of pitch sensor, oriented 
with respect to the m^ frame through a small error angle matrix 


c ± : ci 

Si : Si 


IRA, c 2 = ORA, c 3 = SRA 
IA, g 2 = OA, g 3 = SA 


reference frame fixed relative to gyro case 

reference frame fixed relative to gyro 
gimbal, oriented with respect to the Cj_ 
frame through the gimbal angle matrix 


Angles 




7 g 


'ij 


p . . » 


m. 


■ij 


Euler sequence angles for rotations about the yaw, roll, and pitch 
axes, respectively 

gyro gimbal angle; angle orienting the g ± frame with respect to 
the c ± frame 

elements of the direction cosine matrix orienting the gyro fixed 
reference frame cq with respect to the dq frame 

elements of the direction cosine matrix orienting the gyro fixed 
reference frame ci with respect to the mq frame 

elements of the direction cosine matrix orienting the mq frame 
with respect to the oq frame 


^m/o , ^d/m , l for r iSht-hand rotations in measuring the position of the mq 
() etc. J frame with respect to the oq frame, etc. 


Angular Rates 

j 

to 1 ,o) 2 ,co 3 inertial angular velocity components of the satellite as measured 
in the mq reference frame 

J 

!ww x , angular velocities of the roll, pitch, and yaw inertia wheels, 
respectively, measured in the reference frame 
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Inertias 


1 1 

H 

H- 

C_i. 

1 1 

inertia matrix for the satellite system (with rotating components 
at rest with respect to the satellite) as measured in the mj_ 
reference frame 

1 

H 

H 

O 

O 

1 


0 I 2 0 

principal inertia matrix, that is, inertia matrix in the hj_ frame 

CO 1 
H 

O 

O 


I g 

polar moment of inertia of a single -degree -of -freedom gyro gimbal 

J 1^2 ^3 

vehicle inertia parameters in the b-j_ frame; 

X I 2 - I 3 x II - I 3 x I 2 - II 

1 - T ■ J w 2 - "V ’ d ^ “ T 

ll 12 13 

Angular Momenta 

Hiw 

variable- speed flywheel angular momenta in reference frame 

h i K 

constant -speed flywheel angular momenta in m-j_ reference frame 

Hr 

gyro rotor angular momentum about its spin axis 

Miscellaneous 

c D 

gyro damping coefficient 

M g 

gyro torque generator term 

K 

gyro torque generator feedback constant 

^g 

spurious torques about the gyro gimbal axis 

KcpjKgjK^je 

roll, pitch, and yaw control channel gains 

w 0 

orbital angular velocity 

E* 

torques about the m^ reference frame from sources other than 
wheels or gyros 

T g 

gyro time constant 

Tm 

motor flywheel time constant 
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DESCRIPTION OF THE SATELLITE CONTROL SYSTEM 


A schematic "block diagram of the satellite control system is shown in 
sketch (a). Continuous attitude control of the satellite is through the action 



Sketch (a).- Block diagram of satellite control system. 

of variable -speed inertia wheels nominally alined along three mutually orthog- 
onal, body -fixed axes. These body axes will subsequently be referred to as the 
"marked” axes system and are alined as closely as possible to the principal axes 
of inertia of the satellite. Pitch and roll attitude signals are provided by 
earth horizon sensors. These signals are passed through compensation networks 
before being used as command signals to the pitch and roll inertia-wheel motors . 
The method of measuring angle of yaw makes use of the signal from a single -axis 
gyro oriented so that it senses a component of the orbital angular velocity pro- 
portional to angle of yaw. The output signal from the gyro is passed to compen- 
sation networks before being used as the command to the yaw inertia-wheel motor. 

The part of sketch (a) shown dotted depicts an additional control system 
which operates only intermittently and is used both to despin the inertia wheels 
and to limit the possible attitude deviations of the vehicle. It is not speci- 
fied further since in this report the concern is only with the continuously 
operating inertia-wheel part of the over -all control system. 
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A provision is made in the control system, shown in sketch (a), to inject 
command or "yaw hias n signals into the yaw control channel. In the case of the 
Nimbus satellite, this type of command is used to point a solar array at the sun 
when the sun is not in the orbital plane of the satellite. 


MATHEMATICAL REPRESENTATION OF THE SYSTEM 


In the derivation of the equations of motion of the satellite the following 
assumptions have been made: 

(1) The gravitational field is that which would be produced by an earth 
with a spherically symmetrical mass distribution. 

(2) Coupling between the satellite's translational and rotational degrees 
of freedom is negligible. Thus the rotational equations alone are 
sufficient to determine the angular motion of the satellite relative 
to the orbital path of its center of mass. 

( 3 ) The orbit is circular. 

(4) The torques produced by the gyro rotor and the gyro gimbal are 
negligible. 

(5) The satellite structure behaves like a rigid body. 

The rotational equations of motion, expressed relative to the M marked M 
reference frame, are given below: 


+ ^ 2 ^ 3(133 - I 22 ) + I ! 2 U >2 + Il 3^3 + W!U) 2 I 31 - U>iU> 3 I 2l + l 2 3 (u > 2 2 - U) 3 2 ) 

I 22 (a ) 2 + WiW 3 (lii - I 3 3 ) + 121^1 + 123^3 + W 2 W 3^12 " w l w 2 l 32 + Il 3(^ 3 2 - Ui 2 ) 

J.J- 33^3 + WiW 2 (l 22 - 111 ) + 131^1 + 1 3 2^2 + w 3 w ll 23 - Cx) 3 (j 0 2 I 1 3 + - U> 2 2 ).^ 

Rigid vehicle 


^1W 


W2h 3 w - W 3 h2W 



-J i!i9m/o 


J ili9m/b 

h 2W 

- 

W 3^1W - (^lh 3 w 

+ 3^o 2 < 


-J2l2^m/o 

+ 

J2l20m/b 

_h3W_ 


j^>lh2W - W2hiw_ 


* v 



_ 0 


Variable -speed Gravity torques 

flywheels 


"^haK - w 3^2K' 


' E t i" 

wahiK - wih^ 

4* 


_ w i h 2K “ W2h iK. 


-IX 


v ^ ' 


Constant -speed Torques 

flywheels from sources 

other than 
wheels or 
gyros 


( 1 ) 
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The equation of motion of the combined gimbal and rotor of the single- 
degree -of -freedom gyro is: 


Ig0g + Cjj0g = |^-Ig(C 2 i(i)i + 022^2 + C 23 W 3 ) + H r [(Cn - C3i0g)o)i + ( C 12 ” C320g) u >2 

+ (c 13 - 03302)0)3]^ - Ig[ ( -t^/ m C22 + ^d/m^ 23 ^ 1 

+ ^d/m C21 “ ^d/m 023 ^ 2 + (" 0 d/m C21 + ^d/m 022 ^ 3 ^ 

+ H r [(-t d/m C !2 + 0 d/m Ci 3 )wi + (* d/m Cn - 9 d/m Ci3)w 2 

+ ( _0 d/m c ii + CisTd/m)^ 3 ^ + M g ^ 

where, in general, Mg = -K0g + £Mg. 

The relationship between the inertial angular velocities (wi, 102, and W3) 
and the Euler angles and Euler angle rates are given by the following 
expressions : 

lull = 4k/o cos 0 m/o " V/o sln e m/o cos 9m/ o " ^cA 3111 9m/o sin 0 m/o cos ^m/o 
+ cos 0 m /o Sin trn/o) 

“2 = ^m/o + K/o sin 9m/ o ~ ^o cos 9m/o cos V/o 

W3 = im/o cos e m/o cos 9m/ o + 9m/o sin e m/o _ w o(- cos 0 m/o sin 9m/o cos '•'m/o 
+ sin 0 m/o sin i|r m/0 ) 

Equations (l), (2), and (3) are well known and many derivations of various 
forms of them have been published (see, e.g., ref. 1 ). 

The equations closing the roll, pitch, and yaw control loops are assumed to 
be of the following general form: 


(3) 


Roll h lw - KtpFq/s^gjj./Q 

Pitch h 2 w = KgEg(s)0 S y/ o 


Yaw h 3W = K^F^’ (s) 


(4) 

(5) 

( 6 ) 


where 


G(gyro constant) = 


H, 


K - H r 0) o C32 
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and 


0g(gyro gimbal angle) = 


G[ (Ci 1 s+Ci3O) o )cp m ^ o +Ci2S0 m / o +(Ci3S -Ci 1 0) o )\|r m ^ o 4-Mg-Ci2 w o 1 

1 + T g s 


The development of the gyro gimbal expression from equation (2) is given in 
detail in appendix A. The quantities K<p, Kg, and are the control loop 
gains , and the transfer functions Fcp(s), etc . , represent the combined behavior 
of the horizon sensors, inertia-wheel motors, and compensation networks. It is 
assumed that the transfer functions Fcp(s), etc . , are rational algebraic 
functions of the Laplace transform variable s. 


An important problem that exists for all earth-oriented satellites is caused 
by the set of axes, relative to which the satellite is to be stabilized, rotating 
with respect to inertial space. In the notation of figure 1, the roll and yaw 
axes rotate relative to inertial space with orbital angular velocity, while the 
pitch axis remains fixed relative to inertial space. Therefore, any disturbance 
in the orbital plane which is constant relative to inertial space has oscillatory 
components along both the roll and yaw axes and, clearly, these components are 
90° out of phase. It follows that the rotation of the reference frame itself 
causes coupling between the roll and yaw degrees of freedom. 

In general, the gyro in the control system of sketch (a) will produce a 
signal which consists of a relatively complex mixture of yaw, roll, and pitch 
information, the proportions of each being dictated by the gyro orientation 
within the satellite. It follows that the gyro must tend to introduce coupling 
between the degrees of freedom although, in some cases, it may be sufficiently 
weak to be neglected in the analysis. An important special case of this type of 
control system, and one which will be treated in some detail later in this report, 
is a system in which the gyro is oriented so that its input axis is located in 
the roll -yaw plane of the satellite. In this case it is always possible to a 
first approximation to treat the stability and transient response of the com- 
bined roll and yaw degrees of freedom separately from the pitch degree of freedom. 


Gravity torques play a particularly important role in the dynamics of 
earth -oriented satellites. This can be demonstrated in the following very 
general way. Let G be the vector gravity torque acting on the satellite, and 
co Q the vector angular velocity of the orbital reference frame relative to 
inertial space. The rotational equations of motion of the satellite are: 


GENERAL DYNAMICAL BEHAVIOR OF THE SATELLITE 


Sources of Coupling Between the Degrees of Freedom 


The Roll of Gravity Torques 


H + w 0 X H = G 


( 7 ) 
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It is assumed here that no external torques other than gravity torques act 
on the satellite. The vectors H, G, and co 0 expressed in terms of their com- 
ponents along the o 1; o 2 , o 3 directions (orbital reference frame) are: 


H = O3H1 + 02H2 + O3H3 (8) 

G = Oj_Gi + O2G2 + O3G3 (9) 

CO = -O^COq (lO) 

Thus equation ( 4 ) resolved along the Oi, 03 , 03 directions becomes 

Hi - co 0 H 3 - Gi = 0 (ll) 

H 2 - =0 (12) 

H3 + w qHi ~ G3 = 0 (l3) 


If the gravity torque terms G x and G 3 were missing from the equations, 
then equations (ll) and (13) would represent an undamped sinusoidal oscillation 
of orbital frequency co 0 / 2 jr in the momentum components H x and H 3 . This leads 
to the conclusion that if the satellite system is to be damped, then gravity 
| torques must be present. 


The Effect of Magnitude of the Gain 


The selection of the final compensation networks and gains required for 
j satisfactory dynamic behavior of the class of satellite control systems under 
I discussion, as with most multivariable, multigain control systems, is somewhat 
a matter of trial and error. In this case, however, it is- possible to deduce 
| some important general properties of the system by examining its behavior as the 
gain parameters K<p, Kg, and simultaneously tend to either very small or to 
! very large values. 

i 

j 

When the gains all tend to zero, the system is completely uncontrolled, that 
is, h lv , h 2W , and h 3 -w are all zero. A considerable amount of work has been 
j carried out elsewhere on the subject of the uncontrolled motion of a satellite 
and it is not proposed to repeat this work here except to quote the main results. 

| These are that the energy transferred to a satellite from an outside disturbance 
remains constant in the absence of further disturbances (i.e., the motion is 
undamped) and furthermore the motion is completely governed by the torques 
produced by gravity gradients (gravity torques). 

Perhaps more interesting is the case when the gains IQp, Kg, and are all 
very large. It can be readily seen from equations ( 4 ), (5), and ( 6 ) that the 
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set of variables describing the motion of the satellite body (i.e., 9 m /o> 0 m / o , 
t m / Q ) and their derivatives must all be small compared to the inertia-wheel 
angular momenta h]_v? b 2 v, and h 3 y When the gains tend to infinitely large 
values, the equations of motion (eqs. (l)) become 1 

l23 w o 2 = “^iw + w o^3w + 3 w o 2 J il fPm/b + ^o^sK 

0 = -h 2w + 3w o 2 J2l20m/b (^-5) 

-l2iw 0 2 = -h-3w - w 0 h lw - (l6) 


These equations fall naturally into two independent groups, one containing 
the pitch equation (eq. (15)) and the other containing the roll and yaw equa- 
tions. The pitch equation when integrated yields 


bsw = 3wo 2 J 2 l20 m /bt (17) 

where it is assumed that h2w = 0 when t = 0. This equation indicates that the 
pitch wheel speed increases linearly with time at a rate proportional to the 
pitch control axis misalinement angle 0 m/ q,. The implications of this are 
discussed in a later section on control axes misalinements. 

The roll and yaw equations when integrated yield 

h lw = A cos(w 0 t + e) - h^ + gj 0 I 21 (18) 

h 3 w = A sin(w Q t + e) - h 3K - w 0 (3JiIiqVb " Is3 ' ) ( x 9) 

where A and e are constants depending on the initial conditions. It follows 

from equations (l8) and (19) that there is an undamped sinusoidal oscillation in 

the roll and yaw inertia wheels and they are exactly 90° out of phase with each 
other. This implies that any energy from external disturbances that is imparted 
to the roll and yaw degrees of freedom is transferred entirely to the roll and 
yaw inertia flywheels where it remains indefinitely. Since, with both zero and 
infinite gains, the system has a mode of oscillation with zero damping, then (for 
given compensation networks, inertia -wheel -motor time constant, and gyro time 
constant and provided the system is stable for all finite nonzero values of the 
gains) there must be a set of finite gains which result in the maximum damping 
of the system. 


Equations (l4), (15), and (l6) are derived in the following way: The 

quantities cpm/o* ^m/oj ^m/o and- their derivatives are put equal to zero in 
equations (3) • This yields cj! = 0, = -u) 0 , oo 3 = 0 which are then substituted 

into equations (l). 
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Gyro Orientation Within the Satellite 


Examining the steady-state gyro output signal, 0 ^ steady)' error signal 

used to drive the yaw inertia-wheel motor) will give some indication of the 
influence of input axis orientation on the ability of the gyro to control the 
satellite yaw attitude. An expression for $g( steady) ™ en the vehicle orien- 
tation angles are small can readily be obtained from equation (A 6) by setting the 
Laplace transform variable s equal to zero, that is, 


or 


Hp f i t ZiMp- 

®g( steady) = K - H r w 0 C 32 ^^m/o ' CllU °Vo + hJ 


C 12^0 


0 g( steady) 


-Ciifa/o - C 12 + CisTm/o + (zMg/H r (J 0 ) 
(K/H r w 0 ) - C32 


( 20 ) 


where the quantity AMg represents any spurious torque acting about the gyro 
gimbal axis and the quantity K is the torque generator feedback gain when the 
gyro is operating in the rate mode (K = 0 when the gyro is operating in the 
integrating mode). 

The interesting point about equation (20) is that, provided C 32 is 
nonzero, it shows there is no fundamental difference in the steady-state behavior 
of the integrating and rate gyros. This is due entirely to the fact that the 
frame of reference relative to which the satellite is to be stabilized is itself 
rotating relative to inertial space with orbital angular velocity w 0 . 2 This 
becomes very apparent if oj 0 is made equal to zero in equation (20). The steady- 
state gimbal angle for a rate gyro is then equal to zero or AMg/K, depending on 
whether or not AM g is zero, while for the integrating gyro it is equal to 
(-Ciit m/o - Cis + Ci 3 cp m/o )/-C 3 2 or infinity, depending on whether or not AM g 
is zero. 

An immediate deduction from equation (20) is that, if the gyro signal is to 
be used to control angle of yaw, then Cix must not be zero. An interesting and 
perhaps unexpected implication of this is that a gyro operating in the integrat- 
ing mode, oriented with its input axis along the satellite earth pointing axis 
(m 3 or yaw axis) is unsatisfactory as a yaw control device since Cxj. is zero. 

In addition, in the case of the gyro operating in the integrating mode, it is 
necessary that C 32 no t be zero. 

2 This is also the reason a constant spurious torque acting about the gimbal 
axis of a gyro operating in the integrating mode does not necessarily cause it 
to drift away until it reaches the stops, as would happen if the reference frame 
were not rotating. Instead the gimbal angle reaches the steady-state value 
indicated by equation (20 ) ; that is, if t m /o = 9m/o = ® ’ the va -*- ue 

-C i2 t (AMg/H r td 0 ) 

-C32 
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A desirable, though not strictly necessary, condition is Cj_2 equal to 
zero. A nonzero value results in a steady- state output signal (gimbal angle) 
from the gyro even with tm/o and cfo/o both zero. Although this steady-state 
signal could be removed from the yaw channel by applying a constant bias, it is 
doubtful whether it is worth considering because it is shown in a later section 
that a nonzero C12 has beneficial effects on over-all system stability. 

To meet the zero C12 condition, the gyro input axis must be located in the 
minis or roll -yaw plane of the satellite. In practice it is only possible to 
set C12 equal to zero (or locate the gyro input axis in the did3 plane). 

The value of C12 ma y then differ slightly from zero if the gyro reference 
frame (dq frame) and the control or marked reference frame (mq frame) are mis- 
alined. Thus if C12 is zero, then by equation (A 2 ), is given by 

^d/m Cl1 ' ^d/m 013 ' 

Perhaps the simplest orientation, for the gyro operating in either mode, is 
when C ii is unity, -C 32 is unity, 3 and C 12 and C 13 are both zero. With 
this orientation the gyro input axis lies along the d x or_roll axis of the gyro 
reference frame and the spin axis lies along the negative d 2 axis. The gimbal 
angle expression (eq. (20)) then reduces to 

„ _ - ^m/o “ ^d/m " ^d/irfPm/o + (^g/^r^o) 

8< 6tead y) ‘ (K/ 1 W + 1 l21) 

There may, of course, be objections to this orientation from the point of view of 
over -all system stability. This question is deferred to a later section. It is 
possible, however, to predict a possible objection to this orientation for an 
integrating gyro if it is required to rotate the satellite through large angles 
of yaw, say of the order of 30 ° or more. It can be seen from equation ( 21 ) that 

with K equal to zero the gyro gimbal angle must be approximately equal to the 

angle of yaw, which is clearly incompatible with current gyro gimbal angle limi- 
tations of about 6°. Overcoming this problem, without recourse to special gyros, 
necessitates reducing the value of C n /C 32 from unity to something of the order 
of 0 . 2 . This can be arranged by rotating the gyro input axis in the roll -yaw 
plane so that it makes an angle of about 78° with the roll axis. With this 

orientation, however, the value of C 13 /c 32 is increased from zero to O.98, so 

that the part of the steady- state gyro signal proportional to angle of roll 
becomes about five times greater than that proportional to angle of yaw. Thus, 
not only can the integrating gyro input axis orientation be severely limited by 
the yaw bias requirement, but noise becomes an increasingly important considera- 
tion since the controlling signal in yaw is so much weaker. 

The gyro used in the rate mode does not suffer from problems introduced by 
the yaw bias requirement to the same extent as the integrating gyro. It can be 
seen from equation ( 21 ) that if K is made sufficiently large, the gimbal angle 
can be made small relative to the angle of yaw. 


3 It is shown in appendix A that for a gyro operating in the integrating mode 
it is essential to have C 3 2 negative. 
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ANALYSIS OF THE LINEARIZED EQUATIONS 
Stability of the Linearized Equations 


It is proposed, from the outset, to simplify the subsequent analysis by 
making the following two assumptions: 

(1) There are no misalinements between the marked axes and the principal 
axes of inertia and between the horizon sensor axes and the marked 
axes. 

(2) There is a constant-speed inertia wheel alined along the pitch axis 
only (h^ = h^ = 0) . 

Although the analysis is therefore not completely general, it nevertheless serves 
to indicate many of the main stability characteristics that exist with the type 
of system under consideration in this report. 

The linearized kinematic equations corresponding to equations ( 3 ) are: 


Substituting equations (22) 
izing procedures yields 


= 'Pm/o ' 
= ®m/o 

“3 = i /Q + 

into equations 


u ot m/0 

w o i 


(22) 


“oTm/ o 

(l) and carrying out the usual linear- 


Ii9 m /o - l * w o 2 (l3 - I 2)9 m/0 - Wo(ll - I 2 + I 3)t m/0 + h lw - w 0 h 3W 

" h 2K^m/o + w o c Pm/o) 

1 2 ®m/o + 3 w o 2 (li ■ 4 3 )0 m/o + ^ 2W 

[ l3t m /o + w o 2 ( i 2 - I i)t m /o + «odi - I 2 + I 3)9m/0 + w o h iw + h 3w 

+ ^sK^m/o " w o^m/o^ 




= X) T 2 / (23) 


J 


j The control equations (see eqs. (4), (5)> (6)) become 


hiw - Vcp( s) Vo = 0 (24) 

h 2 w - KgFg(s)0 m / o =0 (25) 
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b-3W + K^(s)[(Ci lS + C l3 u o)9m/o + ^ 12 s ®ui/0 + (^13 s + CiiU 0 )t m/0 


+ AMg - Cl 2 u) 0 ] = o (26) 


where 


*V(s) 


y(s) 

1 + T g s 


(27) 


Equations (23), (24), (25), and (26) form the complete system of linearized 
equations . 


It is convenient, at this stage, to change the time scale from t to r 
where t = 0J o t. This procedure implies that the Laplace transform variable, s, 
in equations (23) to (2 6) must he replaced by co 0 s', where s' is the Laplace 
transform variable corresponding to the new time scale, t. The characteristic 
equation of the complete system expressed in the Laplace transform variable s' 
can then be written in the following form: 


iis' 2 - 4 (i 3 - l 2 ) 


h gK 

w 0 


-di - la + I 3 )s' - 


(ii - I 2 + I 3 )s' + h ^ RS - 

LOQ 


I 2 s' 2 + 3(li - I 3 ) 


I 3 s' 2 + (I 2 - Xi) - if 


0 s' 0 


1 0 s' 


-Kq^cpC^o 3 ' ) 


-KeFgCcOoS’ 


10 0 


1 0 


KyF^(u 0 s' )(cixs' + C13) Ki|rFt|/(w 0 s' )Ci 2 s' K^F^(u 0 s' )(ci 3 s' - Cj.i) 0 0 1 


(28) 


The determinant of equation (28) can be simplified by the following steps: 

(a) multiply row 1 by s' 

(b) add row 3 to row 1 

(c) divide row 3 "by s' 

(d) multiply row 5 by s' 

(e) subtract row 5 from row 2 

(f) multiply column 4 by [(l 2 - Ii) - (bsK/ u o)l 

(g) subtract column 4 from column 3 

(h) subtract row 6 from row 3 
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to 


T „.2 . r,/T t \ . KeFeCwoS' )s' 

I 2 s t + 3(Ii - I 3 ) + 


s’ 


I 1 S ' 2 - 3(I 3 - I 2 ) + II 


s 1 ^ + 1 


II - I2 + I3 + )( C ll sf + C 13 ) 13s' - K^F^ f (oJ 0 S l )Ci 3 S f - Cl! 


-K<pF<p(u)oS ' ) 


0) 0 


(i 2 - ii) - 2* 

w O - 


Expanding the determinant about the second column yields 


= 0 


(29) 


o , KnFn(o) 0 S T ) s T 

I 2 s- 2 + 3di - 13 ) + 0 0 - — = 0 


Ur 


(30) 


and 


I 3 s' - K t F t (w 0 s' )(cl 3 s - C11) 


Iis’ 2 - 3(ls - la) + Ii 


+ 


K<pFcp((jjos' ) (s' 2 + 1) 




(la - Xi) - ^ 
0) o _ 


-4(i 3 - i 2 ) - 


^2K 


(jO r 


(i 3 - i 2 ) + 


h 2K 

U)Q 


s' 2 + ^(woS 1 )(chs + c' 13 ) ( s ' 2 + 1 )| = 0 

(31) 


VJ1 


An examination of equations (30) and (31) reveals that they are independent 
of C 2,2 • It follows that the stability of the system is independent of the angle 
between the gyro input axis and the rKjm 3 or roll -yaw plane. In other words, 
there are no stability improvements to be gained by skewing the gyro input axis 
out of the m.ijn.3 plane. 

Since equation (30) contains only one gain parameter, Kg, the characteris- 
tics of the stability in pitch can be readily obtained by standard root locus 
techniques. On the other hand, equation (31), representing the roll -yaw stabil- 
ity, contains two gain parameters K^p and K^; therefore standard root locus 
techniques do not readily indicate the important features. If, however, the 
additional assumption 4 is made that 

(i 2 - ii) - = 0 (32) 

co 0 

then equation (31) factors into two equations; thus 

I 3 s' - K^Wos'HcIss' - Cii) = 0 (33) 


s 



3ds - la) + Ii 


^ K<pFq)((0 o s' ) (s' 2 + 1) 
O3 0 


(34) 


Equations (33) and (34) each contain only a single gain parameter and loth are 
therefore amenable to treatment by standard root locus techniques . 

While the assumption represented by equation (32) restricts the treatment of 
stability considerably, the class of satellites which satisfy the assumption is 
important. If, for example, the control system does not include a constant speed 
pitch inertia wheel (h 2 j£ = 0), then a satellite configuration satisfying equa- 
tion (32) must have inertial symmetry about its earth -pointing axis. Even if 
equation (32) is not completely satisfied, equations (33) and (3*0 provide a sound 
basis for considering the kind of problems and type of behavior likely to be 
encountered. In fact, it will be shown subsequently that the stability of the 
Nimbus satellite is closely approximated by equations (33) and (3*0 although the 
assumption of equation (32) is not completely satisfied. 

Equations (30), (33), and (3^) can be rewritten in a form more suitable for 
root locus analysis; thus, 


4 The significance of this assumption is that it removes all torques of a 
purely gyroscopic origin from the yaw equation. It is interesting to note that 
this can always be achieved by choosing the appropriate value of h 2K . For 
further insight into the significance of this decoupling parameter, see 
appendix B. 
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x + KgFg(co oS ')s' 
t 0 O I 2 (s' 2 + 3J 2 ) 

K t F t ((o oS ')(Cl3s' - C{i) 

I 3 s' 

IQpFcp(wos' )(s' 2 + 1) 
ojQliS^s ' 2 + 1 + 3J1) 


(35) 

(36) 

(37) 


It is now proposed to examine, in some detail, the implications of equations ( 35)? 
(36), and (37). 


Root loci of equation (35 ) • - First consider the case where no compensation 
is provided; the transfer function F^(oi) 0 s ! ) will then he that of the motor fly- 
wheel in the pitch channel and can he approximated hy a simple lag, that is, 
F^oJqS 1 ) = l/[l + T m (o) 0 s T )]. The type of root locus obtained depends on the sign 
of the gravity torque parameter, J 2 . The root locus diagrams for J 2 positive 
or negative are shown in sketch (b). 5 If J 2 is positive (stable gravity 



Sketch (b).- 


Root locus of 1 + 



1 

1 + T m(wo s ' ) 


(s' 2 + 3J 2 ) 


5 In all root locus diagrams the poles and zeros may he identified hy the 
following code: 

V open loop poles 

M pole associated with motor flywheel 
G pole associated with gyro lag 
G.O zero associated with gyro orientation 

C poles and zeros associated with compensation networks 
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torques), then the pitch channel is stable for all finite positive values of the 
gain parameter Kg and there is a value of the gain which gives the optimum damp- 
ing. When J 2 is negative (unstable gravity torques), the pitch channel is 
unstable for all finite values of the gain. The stability when J 2 is either 
positive or negative can very easily be improved by the inclusion of passive com- 
pensation. If, for example, the compensation is of the often used lead-lag type, 
that is, 

1 + tq(co q s 1 ) 

.1 + a,T e (u) 0 s' ). 


F(a) 0 s’) = 


1 + T m (0J n S' ) 


then with suitable values for t q and cgt q the root locus diagrams can be made 
to appear as in sketch (c). By judicious choice of the lead-lag parameters, the 



Sketch (c).- Root loci of 1 + 


CD 

1 . 

1 + T e (w 0 s' ) 

s 1 

\U3 0 l2/ 

1 + ! ) 

_1 + aT 0 (u) o s' )_J 

(s’ 2 + 3J 2 ) 


0 


optimum damping of the system with J 2 positive can be improved substantially. 
When J 2 is negatave, it can be seen that the stability can be improved so that 
there is a range of values of the gain parameter for which the system will be 
stable even though the damping, in general, will be low. 

Root locus of equation ( 36 ).- The types of root locus diagrams representing 
the solutions of equation ( 36 ) can be divided into two classes corresponding to 
C 13 / 0 and Ci 3 = 0. When C ' 13 / 0, the equation can be written in the form 
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1 


(38) 


(Mk) v . 0 ,.) t- - _ 


= 0 


and when 0^ = 0, in the form 


1 + 


t C ll\ Pf(wo s ' ) 


Is 


= 0 


(39) 


First consider the case where no compensation is provided. The transfer 
function F^(a) 0 s , ) can he written in the form 


tyMos’ ) = 


1 + T m( w o s ') 1 + T g( w o s ') 


where T m is the motor flywheel time constant and Tg is the gyro time con- 
stant. The root loci of equation ( 38 ) are shown in sketch (d) for Cii/Cia 





(4) 


j" 

~C,3 



~Ci7 T m"o 



GO. M G 


y 



1 a 


Large r g 


Sketch (d).- Root loci of 1 



1 1 

1 + T m( w o s ' ) 1 + T g (w 0 s’ ) 


[s’ 


(CliMs)] 

s' 


= 0 
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negative and positive. The case where C11/C13 is negative is subdivided into 
two cases corresponding to -C'^/Cxx < u) 0 T m and -Cis/Cix > co 0 T m . The root loci 
of equation (39) (i.e., Ci 3 = 0) are shown in sketch (e). Two values of Tg are 



considered: one very small corresponding to a rate gyro, the other very large 

corresponding to an integrating gyro. The small value of r g is less than both 
u) 0 T m an( 3- IC13/C11I and the large value is greater than both u) 0 T m an( ^ I C13/C11 I • 

The root loci diagrams shown in sketches (d) and (e) reveal the important 
fact that the best damped systems are those with small gyro time constants 
(i.e., rate gyro) and with gyro orientations such that Cii/Ci 3 is negative. 

To obtain a geometrical interpretation of the latter of the two requirements, 
consider the case where the Cj_ frame is oriented to the m^ frame by a stand- 
ard pitch, yaw, roll Euler angle sequence (similar to that used in appendix D to 
orient the c± frame with respect to the d^ frame). If the pitch, yaw, and 
roll angles are designated 0 c / m , ^c/rru Tc/rrp then 


CI3 = sin e c /m cos »c/m 
°ii cos e c/m cos * c / m 


(40) 


Thus a negative value of C11/C13 means that 0 c / m must lie in either the 
second or fourth quadrants. It makes no difference which of these two quadrants 
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is selected, provided a K\|/ of the correct sign is used. Sketch (f) gives a 
pictorial idea of a gyro oriented with its input axis in the fourth quadrant, 
when Kxjf must he positive for stability. 



To obtain a gyro signal which has the component proportional to the angle 
of yaw a maximum, it has been shown that the input axis must be located 

along the dx axis (see eq. (21) and related discussion). The advantage of 
obtaining the strongest yaw signal is, of course, that the influence of noise 
and signal threshold limits is minimized. With an orientation of this type, C u 
is unity and C 12 and C 13 are both zero. The expressions for Cx 3 and Cxi 
obtained from equation (A2) are 

C13 = ~®cL/m 
Cii = 1 

Therefore 


p . " d/m 
^11 

If 0q/ m is less than or equal to zero, then Cx 3 /Cxi is greater than or equal 
to zero, and it can be seen from sketches (d) and (e) that the system has an 
oscillatory instability for values of the gain parameter greater than some crit- 
ical value. To guard against the possibility of an instability due to gyro mis- 
alinement (i.e., %/ m ) ? the gyro input axis must be oriented at an angle, 9 c / m , 
as in sketch (f), so that 


tan 0 c / m < 0 
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or, since tan 0 c / m need be only slightly less than zero, 

®c/m < 0 

Since 0 c / m = 0c/d - ^d/m* 6 "this criterion for stability may be written as 

®c/d ^ ®d/m 

The least value of 0 d/m may be written as -| 0 d/m | ^ where j 0 d/m | ^ is 
the estimated maximum numerical value for the misalinement error angle 0 d / m . 
Hence the criterion for stability becomes 

®c/d ^ "| ®d/m|max (41) 


In the case of the integrating gyro (large T g ), it can be seen from 
sketch (d) that, in general, a large value of -C X3 /C XX yields better stability 
characteristics than a small value. In other words a large value of -0 c / m is 
desirable (direction of -0 c / m as in sketch (f)). It is interesting that this 
condition is compatible with the integrating gyro orientation conditions (estab- 
lished earlier in the section dealing with gyro orientation) necessary to meet 
the yaw bias requirement. 

Even if the gyro orientation is favorable and the time constant small, the 
damping factors of the oscillatory roots at high gain may be unacceptably low 
for some applications. This can be improved by adding lead-lag compensation in 
a manner similar to that indicated in sketch (c) for the pitch channel. 

Root loci of equation (37) - - Again consider first the case where no 
compensation is provided. The transfer function Fq>(u 0 s’ ) is then the simple 
lag function assumed to represent the motor flywheel dynamics. The type of root 
locus obtained depends on the sign of the gravity torque parameter J x . The root 
locus diagrams for J x positive and negative are shown in sketch (g) . The 
system with J x positive ha s a conve rgent mode of oscillation lying in the 
frequency band co 0 /2n: to gj 0 n/ 1 + 3 J 1 / 2rt ; whereas with J x negative, such that 
0 > J x > -l/3 , the mode of oscillation in this frequency band is divergent. If 
J x is in the range -l/3 > Ji > -1, the divergent mode may be oscillatory or 
nonoscillatory depending on the value of K^. If J x is zero (no gravity 
torques), the poles and zeros on the imaginary axis become coincident and there 
is an undamped sinusoidal oscillation of orbital frequency u 0 /2n: for all values 

of the gain K<p. This is in agreement with the general conclusions given earlier 
in the report. It should also be noted, again in agreement with earlier con- 
clusions, that, as the gain is made very large, the system tends to possess an 
undamped sinusoidal oscillation of orbital frequency. 


6 The negative sign for 0 d / m occurs because 0 c /m an ^ ®c/d are rotations 
about the c 2 axis while 0 d / m is a rotation about the d 2 axis and by the 
convention adopted in appendix D, c 2 = -d 2 . 
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It is shown in appendix B that equation (3*0 governs the transport of 
roll -yaw angular momentum between the satellite and its orbit. (When the trans- 
port of angular momentum is from the satellite to the orbit, the process is 
usually termed "momentum dumping.") It is clearly desirable to maintain the 
angular momentum of the satellite (including the inertia wheels) as small as 
possible within the framework of other satellite performance requirements. It 
can be seen from sketch (g) that the best system from the momentum dumping point 
of view is the one with the largest positive value of J x (e.g., a long slender 
cylinder with its axis earth pointing). With J x negative, momentum dumping is 
not possible with an uncompensated inertia-wheel control system. It is not 
possible to materially improve a system with J x between 0 and -l/3 even with 
very general forms of passive compensation. It is theoretically possible to 
improve a system with J x between -l/3 and -1 by cancelling the unstable pole 
with a zero, but it is doubtful that the resulting system would be practical. 
Even with J x positive it is very difficult to materially improve the maximum 
damping (given by point A of sketch (g)). Furthermore, the gain corresponding 
to the maximum damping is relatively small and for many systems may not be com- 
patible with the demands of other performance requirements. If these require- 
ments dictate a high value for Kp then the momentum dumping capability will be 
extremely low and there is little that can be done to improve it. However, the 
fact that one of the roots of the characteristic equation is almost undamped 
does not necessarily mean that the system is unsatisfactory. To demonstrate 
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this point, consider such a system at a time sufficiently removed from the time 
of a disturbance that all the transients, with the exception of that associated 
with the lightly damped mode, have died away. The system will then he oscillat- 
ing at close to orbital frequency with approximately constant amplitude, so that 
the Laplace transform variable s can be set equal to jto 0 (j = >7-1) in the equa- 
tions governing the system. Thus the roll control equation (eq. (24)) becomes 


h rw _ 


c Pm/c 


= Kq>Fq)(ju) 0 ) 


(42) 


or equating the modulus on both sides of equation (42) 


h 


1W 


Vo 


Ih 


iw i 


4Wo I 


- Kepi Fq>( jw 0 ) 


(43) 


We now define | h lw | may to be the maximum inertia-wheel angular momentum, at 
which the coarse control system is activated to desaturate the inertia wheel, 

| Tm/ol max spec and | Vo I max spec to be the specified deviations 

of roll angle and roll rate, respectively. If Tm/omax and ^m/oimax are 
the maximum values of cp m/Q 9m/ o corresponding to | h^lmax^ then the 

following are specified: 


Vo 

max < 

Vo| 

4Wo 

max ^ 

^m/o 


max spec 
max spec 

hut since the motion is sinusoidal with circular frequency u) Q 


(44) 

(45) 


'Vo max w o c P m /o max 
Thus inequalities (44) and (45) may be restated in the following form: 

I TWoJmax spec 


9, 


m/o 


max 


< -^greatest of ] 9 n - /0 


max spec 


and 


tor 


or from equation (43) 


(46) 


^greatest of • 

l tLrwImax 

OYl H 

“ol^ivlmax l 

^cp( jto 0 )| 

4*m/o 

CU.1LL 

max spec 

Fcp( ju 0 )| 

Vo 

max spec 


(47) 
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Therefore if is large enough to produce a system with small damping in the 

sense discussed previously, then inequality (4-7) is a necessary condition that 
the system will he satisfactory in meeting the specifications. In a similar 
manner to the above, an inequality connecting with | h 3 y[ m a y . ^m/o jmax spec 

and i tm/n | may spec can derived which must also be satisfied. It is not pro- 
posed to go into the details of this since the expressions involved are long and 
not particularly revealing. 

If, in order to meet the specifications, a very high value of Kp is 
required, then it may be desirable to improve the damping factor of the pair of 
oscillatory roots located on the branch marked B of the root locus diagram shown 
in sketch (g). This can be accomplished by the addition of simple lead-lag com- 
pensation applied in a manner similar to that indicated in sketch (c) for the 
pitch channel . 

The effect on stability of (l 2 - lx) and h 2 K - - If the system does not 
satisfy equation (32), even approximately, then it is necessary to consider the 
complete characteristic equation (eq. (31)) in determining the stability of the 
system. Unfortunately as with most multivariable problems of this kind it is 
very difficult, if not impossible, to deduce any general properties of the sys- 
tem by methods currently available. Several specific examples have been evalu- 
ated in which it appears that small negative values of (l 2 - Ii)/l3 = J3 an< l 
small positive values of h 2 K/u) 0 I 3 result in a system with a maximum momentum 
dumping capability which is slightly better than that which would be predicted 
by equation (37), although always at the expense of the damping of other modes. 
Larger values of -J 3 and h 2 K/o) 0 I 3 eventually result in a decrease of the 
momentum dumping capability below that predicted by equation (37)* With very 
large values of h 2 K/oj 0 I 3 , the dynamics of the constant speed wheel tend to 
dominate the over-all stability of the system and always result in an almost 
undamped oscillation of orbital frequency which tends to create large amplitude 
roll and yaw deviations. Thus for the type of satellite system considered here 
any improvement due to the use of a constant speed wheel is always small and 
probably not worth the additional mechanical complications. 

A possible way of improving the roots of equation (37) for the cases where 
J 3 is nonzero, but small, is to replace the imaginary poles by the correspond- 
ing true open loop poles derived from the open loop characteristic equation 

s' 4 + s ,2 (l + 3Ji + J 1 J 3 ) + 4JiJ 3 = 0 (48) 


The Influence of Axes Misalinement 

Two types of axes misalinement can occur and can influence the behavior of 
the satellite significantly: 

(a) Angular misalinements of the axes system selected as the reference for 
control information (known as the marked or control axes). 

(b) Angular misalinements of the sensors relative to the marked axes system. 
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It is usual to attempt to aline the control axes so that they are coincident 
•with the principal axes of inertia of the satellite. If some orientation other 
than this is selected, then in its equilibrium attitude, the satellite generally 
has torques acting on it as a result of the action of the gravity gradient. 

These gravity torques have to he countered by the control system and the only 
way this can be achieved is to continuously change the speed of the inertia 
wheels. This situation is undesirable for two basic reasons. First, there is 
a tendency to reduce the life of the inertia-wheel units and, second, an addi- 
tional burden is placed on the coarse control system (and its energy supply) 
through the increased likelihood of the inertia wheels attaining saturation speed. 

The difficulty of accurately determining the location of the principal axes 
of inertia of the satellite makes it almost impossible to completely avoid the 
problems outlined above. In addition, the uncertainty of the location of the 
principal axes can introduce important stability problems, as will be demon- 
strated subsequently. Opinions differ somewhat on the question of attainable 
accuracy with which the principal axes can be determined, and, certainly, the 
type of configuration is an important factor. It appears, however, that the con- 
trol system design should allow for the possibility of a control axis misaline- 
ment of at least 3° and possibly as much as 5°. 

Misalinement of the sensors can result from manufacturing errors and 
structural distortion. Their effect is to produce both attitude errors and 
inertia-wheel speed biases. In general, however, the influence of these mis- 
alinements is largely overshadowed by that due to the control axes misalinements . 

It is possible to obtain a quantitative idea of the influence of misaline- 
ments by considering the steady-state equilibrium of the system. The steady- 
state equations governing the static equilibrium of the system can be derived 
from equations (l), (3), (4), (5), and (6). These equations take on a simplified 
and more meaningful form when linearized, a process which is justifiable in the 
present circumstances since no external torques are assumed and the misalinement 
angles are, in general, small. The equations are given below. 


Roll 




,s 

^ 3 W 

I1W0 

( 49 ) 

Pitch 





= J 20m/b (50) 

Yaw 



( ^2K , T ^ , h iw 

\ W 0 I 3 ' V m/0 13^0 


= J 3V/b (51) 
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control 

o' 

o 


v,S 
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w 0 


-IQpcPsx/m 
100 ( 52 ) 

Pitch 

control 


%^m/o 

- h|w 


( 53 ) 

Yaw 

control 

K\|/Ci 3 cp^/ 0 



h|w 

w o 

= 0 (54) 
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where the suffix s means the steady-state value of the variable. It should be 
noted here that misalinements of the gyro orientation enter the equations only 
through the terms Cix and Ci 3 . Thus gyro misalinements do not of themselves 
introduce any steady-state attitude angles and inertia-wheel speeds although they 
do influence the value of those derived from other sources of misalinement . 

It follows directly from equations (50) and (51) that in the case of pitch 


4/0 = e m/b <55) 

and 

h 2 w = K e(0m/b + e sx/m) (56) 

provided 

j 2 t 0 


Equation (55) merely states that, in equilibrium, the principal axes of inertia 
in roll and yaw always point along the roll and yaw axes of the orbital reference 
frame (i.e., 1d 1 = o^ and b 3 = o 3 ). As the pitch gain Kq is increased, equa- 
tion (56) indicates that the steady-state -pitch wheel speed increases. If the 
gain is high, however, the wheel speed will probably saturate before it reaches 
the steady-state value given by equation (56) and will therefore have to be 
desaturated by the coarse control system. For high values of the gain the rate 
of increase of wheel speed with time tends to the asymptotic value (corresponding 
to an infinite gain) given by equation (17). Thus for high gain systems, equa- 
tion (17) provides a rapid approximate estimate for the time interval between 
wheel desaturations for a given pitch misalinement angle. 


The steady- state roll and yaw angles (cp^/ 0 and ^m/o) an(i associated 
inertia -wheel angular momenta hf w and h 3w can be obtained by solving the simul- 
taneous equations (49) > (51), (52), and (54). The general results are not pre- 
sented here since the expressions are not particularly revealing. When the gains 
K\|f and Kcp are infinite, the general expressions reduce to the following 


g 

c Pm/o = _c Psx/in 
hiw = I3J3W0 ( ^sx/m 5^ + 

^m/o - C ^ sx / m 

^3w = -Ii w o^i( c Psx/m + ^m/b^ + ^sK^sx/m 




m/b 


- h 2K5P 


ci 

sx/m c'. 


13 


11 




(57) 
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The important properties indicated hy equations (57) are that the inertia-wheel 
speeds are "bounded (note that the pitch wheel speed is unbounded when K@ tends 
to infinity) and the steady-state angles cp^/ Q and t^ /0 are dependent only on 
the sensor misalinements which are usually very small (of the order of l/lO°). 
These facts illustrate a general point noted earlier in connection with stability, 
namely, that misalinements in pitch have' a much more significant effect on the 
behavior of the system than have misalinements in roll and yaw. 


The Influence of Control System Failure 


The desired lifetime of a satellite is usually of the order of months and 
in some extreme cases, years. The multiplication of components to achieve the 
desired reliability is often very costly in weight. In some cases, however, it 
is possible, after a component fails, to rearrange the operation of the remaining 
components so that stability is maintained with some degradation of performance. 
Clearly, it is desirable to exploit such circumstances to the maximum since this 
may provide a way of improving system reliability for only a small weight penalty. 

The stability of the roll and yaw degrees of freedom can always be arranged 
to be dependent on each other. The condition for this is that the system must 
not satisfy equation ( 32 ), or the system parameter must satisfy the equation 


- ~ ^ 1 ^ = K (where K £ 0) 

I 3 Cx) 0 l3 

This raises the possibility, in the event of a roll or yaw channel failure, of 
activating a constant speed wheel to adjust the value of K so that the result- 
ing system is stable. Furthermore, if the value of h 2 £ required to produce the 

desired value of K is not too large, there is the possibility of using the 

pitch control inertia wheel to provide the value of h using the remainder of 

its speed capability for pitch control in the usual way. It should be noted that 

the coupling between any of the degrees of freedom can be changed by changing the 
values of h^, h 2 K> and h 3 £. For example, it may be possible to minimize the 
effect of a pitch channel failure by selecting appropriate values for either hjja 
or h 3 K or both, which can be provided by the roll and yaw control wheels, respec- 
tively. However, in all but the case considered here (i.e., varying h^) th e 
analysis appears to be rather intractable from the point of view of yielding 
general conclusions and the feasibility of these schemes can only be determined 
for specific satellite configurations. 

In the subsequent discussion on the influence of a control channel failure, 
it is assumed that the system prior to failure has those design features which 
have been shown to be desirable -in the previous discussion. In particular the 
following are assumed: 


28 


J x > 0 


^<0 

Cis 




Trr small 

G> > 


(58) 


Failure of the yaw control channel.- The stability of the roll -yaw motion 
with the yaw control channel inoperative may he obtained from equation ( 31 ) by 
setting equal to zero. The resulting equation may be written in the follow- 

ing form, suitable for a root locus investigation: 
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(59) 


It is shown in appendix C that the best open loop stability (i.e., K<p = 0) is 
obtained when h 2 ^;/u 0 satisfies the condition 


[both (l 2 - lx) and 4(l 2 - I 3 )] 

w 0 “ 


The open loop poles are then all located on the imaginary axis and are, of 
course, symmetrically disposed relative to the real axis. The two poles on the 
positive imaginary axis are located so that one is always equal to or less than 
j and the other always greater than or equal to j \TC , where C is given by 
the expression 


C 
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Sketch (h).~ Root loci of equation (59)- 


The three types of root loci that can he obtained are shown in sketch (h) where 
F C p((jJ 0 s T ) has been assumed to be a simple lag function approximating the motor 
flywheel dynamics. It can be seen that a stable system is achieved only if there 
is an open loop pole on either side of each of the zeros located at the points 
s T = ±j (i.e., sketch (h)l and 2). This can be arranged by selecting an appro- 
priate value for h 2 ^/a) 0 (for example, any value for which C > 1 would be satis- 
factory) . The value of K^/I^q corresponding to maximum damping of the system 
represented by either 1 or 2 of sketch (h) is always relatively low (i.e., com- 
parable with its value at point A of sketch (g)). Hence if the original system 
has high gain, this method of overcoming the effect of a yaw channel failure is 
not very satisfactory, since two almost neutrally stable roots occur, one at 
s T = ±j, the other at s T = 0. Whereas the neutrally stable root of orbital 
frequency (s f = ±j) can be tolerated in the original system (since most of the 
system angular momentum is always stored in the inertia wheels) after the system 
has a failure there is no yaw wheel to store the yaw angular momentum which must 
therefore appear in the form of large yaw angles and yaw rates. 


Failure of the roll control channel . - The characteristic equation for a 
system with a roll channel failure is given below (eq. ( 31 ) with K^p = 0) 
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The zeros of equation (60) may be obtained in the form of a root locus diagram 
with a gain parameter -[J 3 - (h^/wola)] I 3 /I 1 Cxi/Ci 3 . Thus the root locus 
equation for the zeros is 
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A typical pole -zero configuration corresponding to equation (6l) is shown in 
sketch (i), where conditions ( 58 ) have heen assumed to apply. 

Since F^(u) 0 s T ) consists only of 
poles in the left half plane and the 
most stable open loop poles of equa- 
tion ( 60 ) are all on the imaginary axis, 
it follows that the existence of a zero 
on the positive real axis (zero A of 
sketch (i)) implies the system repre- 
sented by equation ( 60 ) is always 
unstable. It must, therefore, be con- 
cluded that the scheme as outlined is 
not, in general, satisfactory in the 
event of a roll wheel failure. It is 
possible to visualize a modification to 
the scheme whereby additional compensa- 
tion is switched in after a roll wheel 
failure. Such modifications are, how- 
ever, outside the scope of this report. 


A 

-e- 
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Sketch (i).- Typical pole-zero con- 
figuration for equation (6l). 


Steady-State Behavior at Large Angles of Yaw 


For some purposes it may be required to rotate the satellite to some yawed 
position while maintaining the pitch and roll attitudes close to zero. One 
example of this is the previously mentioned yaw bias requirement for Nimbus. Yet 
another example could be a requirement to point, in a given direction, an antenna 
which is rigidly fixed to the main body of the satellite. An important property 
of the type of system under discussion is that if the yaw angle is large, then 
the nonlinearities and couplings inherent in the system produce steady-state 
pitch and roll attitudes which differ from zero. (it is assumed here that no 
provision is. made to inject command signals to the roll and pitch channels.) The 
error involved in these angles is dependent on the angle of yaw, the inertia dis- 
tribution of the satellite and the roll and pitch channel gains. In order to 
calculate the steady-state roll and pitch angles, it will be assumed that they 
are sufficiently small that their products may be neglected. The steady-state 
inertial angular velocities from equations ( 3 ) are 
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(62) 


“I s = " w o sin \|f® /o 

co 2 s = _ Uo cos I 

W3 s = -(o 0 [-cp^ /o cos l£ /0 + 0^ /o sin 4 /0 ] 

where the superscript s means the steady-state value of the variable. Sub- 
stituting for W! S , W2 S j and ^3 S from equations (62) into equations (l) yields 

(Oo 2 COS ^ /o [-cpm/o cos *m/o + 0 m/o sin ^ /o ](l 3 - la) 

= co 0 cos t^/ohlv-ojoC-cPm/o cos *m/o + 0 m/o sin +5/ 0 KhIw + h sK ) - 3^ 0 2 (l2 - I 3 )cp ® /o 

(63) 


sin Hft /o [-q£ /o cos fS /Q + 0 s /o sin ^s /o](l] 


I 3 ) 


= w q( _c Pm/o cos ^m 


(64) 

w 0 sin cos tm/o(-^-2 - Ii) = w 0 (h2w + h2^)sin - co Q hi W cos tm/o 

(65) 


where misalinement errors have been assumed to be zero. The result of multiply- 
ing equation ( 63 ) by sin equation (64) by cos \|r® /0 , and adding is, after 

some rearranging, 

[ _cp m/o cos ^m/o + 0 m/o sin Co^'^sin ^ /Q cos t° /0 (l 2 - I x ) 

+ w Q (hl w + h 2K )sin \|r® /0 - to Q hi W cos i|r® /0 ] 

= -3w 0 £ [sin t® /0 (l 2 - I 3 )q £ /0 + cos l£ /0 (li - I 3 )0® /o ] (66) 

hence from equation ( 65 ) it follows that equation (66) reduces to 

sin ^m/o( l2 - I 3)9m/o + cos ^m/ o ( 1 1 " = 0 ( 6 7) 

The steady-state roll and pitch inertia-wheel angular momenta hx W and h! w are 
given by equations (4) and (5), that is. 
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h ?w = v4/o 


( 68 ) 


S g 

h 2w = K00 m/O 


(69) 


where, again, misalinements are assumed to he zero. Using equations (68) 
and (69), equation (65) can he written in the following form: 
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(70) 

Equations (67) and (70) form two simultaneous linear equations in 9^/o an<3 - *?m/o 
whose solution is 
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sin 2 tm/ 0 ( I 2 - 13)% + cos 2 tm/ 0 ( I i - *3)^ 


(72) 


An important property illustrated hy equations (71) and (j2) is that if 


(I 2 - I 1 )cos t®/o - ^ = 0 (73) 

then 0^/0 an< i 9m /o are zero independent of the values of the gains K@ and K<p. 
For some mission requirements this could he a very valuable property and could 
he achieved for a given yaw angle *I/o simply hy selecting the appropriate 
value of h 2 K* (Note that equation (73) is identical to equation (32) when 
*&/o = 00 Another way of reducing 0^/0 an & Vm/o i ncrease the system 

gains Kq and K^p provided this does not violate any other requirement or lead 
to an unsatisfactory dynamic behavior. 
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COMPUTED RESULTS FOR A PRACTICAL EARTH -ORIENTED SATELLITE 


Satellite inertias and "basic control system data are given in appendix D. 
These data were typical of the Nimbus meteorological satellite at one stage in 
its development. The inertia distribution is favorable (Ji >0 and J 2 > 0 ) in the 
sense that the roll and pitch gravity torques are stabilizing. Identical lead- 
lag compensation networks are used in the pitch and roll channels . No compensa- 
tion is used in the yaw channel, and no constant speed inertia wheels are used. 

Studies were carried out using both digital and analog computers. Calcula- 
tions on the digital computer were confined to investigations into the stability 
of the linearized equations describing the system. In the analog simulation, use 
was made of the complete equations with all the coupling and nonlinearities 
included. The analog study was used to determine the transient response of the 
system, particularly following a yaw command, and to check the validity of the 
stability derived from the linearized equations. 

It has been demonstrated previously that the smaller the gyro time constant, 
the better the stability characteristics of the system and, furthermore, that it 
is necessary to incline the input axis of the gyro relative to the roll axis in 
order to avoid the possibility of instabilities being caused by gyro misaline- 
ments. It will be demonstrated subsequently that additional gyro inclination to 
the roll axis is necessary to avoid instabilities which can be caused by mis- 
alinements of the "marked" axes relative to the principal axes of inertia. In 
order to meet these desirable features, the basic control system studied is 
assumed to have a gyro with a time constant of 0.2 sec oriented with its input 
axis at an angle of -6° to the roll axis, that is, cp c /^ = 0°, = -6°> and 

\|r = 0 ° (see appendix D) . 


Influence of Control System Gain 


Root locus diagrams for the roll and yaw degrees of freedom of the basic 
system are shown in figures 2 and 3 > respectively. It has been assumed that 
equation (32) is approximately satisfied, which, in this case, amounts to the 
statement that the pitch and roll inertias are approximately equal. The advan- 
tage of using an approximate solution of this kind for the roots of the roll -yaw 
characteristic equation is that it enables the stability of the roll and yaw 
degrees of freedom to be examined separately, thereby greatly facilitating an 
understanding of the influence of various system parameters, in particular, the 
system gains. The validity of the approximation in this case can be seen from 
table I, which shows a comparison of the approximate roots obtained from equa- 
tions (36) and (37), an improved approximation obtained by using the true open 
loop poles on the imaginary axis of the roll diagram, and the exact roots. The 
comparison is given for two sets of the gains K<p and K\|r. The set with the 
smallest values (i.e., Kp = 0.3 ft -lb -sec/radian, Ity = 300 ft-lb-sec 2 /radian) 
results in a system with the greatest damping and one which therefore makes the 
greatest use of gravity torques to remove angular momentum from the vehicle (see 
fig. 2 , also sketch (g) and related discussion). This system will subsequently 


be referred to as the "low gain system. " The other set of gains given in 
table I results in a system which, for all practical purposes, is undamped (see 
fig. 2 ) . The gains are, however, sufficiently high to cause most of any initial 
angular momentum of the system to pass quickly to the inertia wheels. As a 
result the vehicle body angles and rates, in the residual undamped oscillation, 
are very small. This system will subsequently be referred to as the "high gain 
system." It can be seen from table I that the results of both methods are prob- 
ably accurate enough for most engineering purposes. The second method of approx- 
imation produces better results for the low gain system while the first method 
produces better results for the high gain system. It is interesting to note that 
for the low gain system the coupling between roll and yaw introduced by the 
unequal roll and pitch inertias changes a real roll root and real yaw root into 
a complex pair of roots (see table I). 

To complete the stability investigation, a root locus diagram for the pitch 
degree of freedom is shown in figure 4 . This is similar to that shown in 
sketch (c) and discussed earlier. In order to maintain a balance in the stiff- 
ness of the control channels, the pitch gain Kg for the low gain system is 
made equal to 0-3 ft-lb-sec/radian and that for the high gain system 
76.4 ft-lb-sec/radian. It should be noted that, for both low and high gain 
systems, the channel gains have been arranged so that one unit of satellite body 
angle about any degree of freedom produces, in the steady state, approximately 
the same value of the angular momentum of the corresponding inertia wheel (i.e., 
Kg = K<p = to 0 K\j;) • 


A comparison between the transient motion of the high and low gain systems 
when released from two different initial conditions is given in figures 5 and 6. 
It is immediately apparent that the transients of the low gain system are com- 
pletely dominated by the lightly damped (£ = 0.13) mode of approximately orbital 
frequency. Satellite body angle, body rates, and inertia- wheel speeds all damp 
to zero although it takes several orbits for this to occur. Perhaps one of the 
most striking features of the low gain system behavior is the absence of any 
sharp peaks in any of the variables describing the system response. In contrast 
the high gain system, especially for initial conditions involving satellite body 
angles, exhibits rapid initial transients with relatively high peak values of 
body rates and inertia-wheel speeds - a behavior fully in accord with the charac- 
teristic roots of the high gain system (table I). When the rapid transients have 
died away, there remains, in the roll and yaw degrees of freedom, the residual 
undamped low-frequency mode. This shows up particularly in the inertia-wheel 
speeds. It does not show on the traces of the satellite body angles and body 
rates because the gain of the system is so high that the amplitude of the mode 
is below the resolution of the recorder. 

It can be seen from figures 5 and 6 that there is only light coupling 
between both pitch and roll and pitch and yaw. To illustrate this further and to 
indicate the type of pitch response obtained, figure 7 shows the transients 
following a release of the system from an initial condition of 5 ° of pitch. As 
might be expected, the response of the low gain system is much slower than that 
of the high gain system although the damping remains very good. 

There is no doubt that the high gain system provides a measure of control 
over the body rates and body angles which is far superior to that provided by the 
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low gain system. The principal objection to the high gain system is that accumu- 
lated system angular momentum cannot he hied off through the action of gravity 
torques. Instead, the only way of removing angular momentum from the system is 
hy the expenditure of stored energy through the action of the coarse control 
system. 

There remains the possibility that some control system with intermediate 
gains may provide some of the good response characteristics of the high gain 
system while maintaining sufficient damping of the low frequency mode to enable 
some degree of momentum dumping capability. An example of such an intermediate 
gain system is shown in figure 8. Here the gains have been set at 
YLq = K<p = 1 ft -lb -sec/radian and K\|r = 1000 ft-lb-sec 2 /radian. The damping factor 
of the low-frequency mode is reduced from 0.13 for the low gain system to 0.055* 
Thus the low-frequency mode takes over twice as long to damp to a given fraction 
of its initial value; however, the response of the system is certainly improved. 

Whether or not a compromise such as this can be regarded as satisfactory 
depends entirely on the design specifications. If the body angles and body 
angular rate specifications are not too stringent, then such a system could well 
provide an acceptable solution and one which requires the minimum amount of 
stored energy within the satellite. 


Influence of Gyro Time Constant and Orientation 


It has been suggested previously that, in general, it is desirable to use 
a gyro with a small time constant oriented to give a small positive value of 
-Cis/Cii.. The basis for the argument is that these conditions result in the 
most stable roots of the yaw characteristic equation (eq. (36)). It has also 
been indicated that a system using an integrating gyro (very large time constant) 
can be made to yield reasonably stable roots of the yaw characteristic equation, 
provided the gyro input axis is oriented to give a large positive value of 
“C13/C11. In general these roots are not as stable as those for the system using 
a rate gyro. However, if these roots are such that they do not have an over- 
whelming influence on the transient behavior of the system, then it may be pos- 
sible, especially with low control channel gains, to devise a system, using an 
integrating gyro, with a performance comparable to that using a rate gyro. The 
reason for believing this may be possible for the low gain system is that the 
transient response of this system has been shown to be primarily dependent on the 
low frequency oscillatory mode associated with stability in roll. 

To check whether or not the above conjecture is valid, the system was 
changed in two stages. First the rate gyro (time constant of 0.2 sec) was 
replaced by an integrating gyro (time constant of 2000 sec), gyro orientation 
remaining unaltered at (0°, -6°, 0°). The main observable effect on the tran- 
sient behavior of the system is an apparent increase in the coupling between roll 
and yaw in the sense that a disturbance introduced into one degree of freedom 
produces greater amplitude transients in the other (cf. figs. 5* 6, 10, 11, 

12, 13, Ik). The initial nonlinear behavior displayed in most of the traces 
involving the high gain system are due to inertia -wheel motors attaining maximum 
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speed or maximum torque. The explanation for the coupling is that the roots of 
the yaw characteristic equation, which in the basic system are very damped, have 
had their damping reduced to the point where they materially influence the tran- 
sient behavior of the system. The second change to the system involved a reori- 
entation of the integrating gyro so that its input axis was inclined at an angle 
of -70° to the roll axis. Thus the value of -Cis/Ci! is increased from 0.105 
to 2.748 and changes the yaw root locus diagram from the type shown in 
sketch (d)4 to that shown in sketch (d)3» The effect of this change on the tran- 
sient behavior of the low gain system is shown in figures 9 and 10. It can be 
seen that the improvement in the damping of the roots of the yaw characteristic 
equation correspondingly improves the transient response of the system. In fact, 
a comparison of figures 9? 10 and 5? 6 shows that there is little to choose 
between the transient behavior of this system and that of the basic system 
employing th§ rate gyro. The effect of the system change on the high gain 
system is shown in figures 11 and 14. Although the change of gyro orientation 
greatly improves the transient behavior of the system, it is still considerably 
worse than that of the basic system. This is particularly evident when the 
initial conditions are in the form of angular rates, where a relatively poorly 
damped oscillation and a very slow convergence tend to dominate the motion. 

To demonstrate the influence of gyro orientation on the behavior of the 
system using a rate gyro, the orientation was changed from (0°, -6°, 0°) to 
(0°, 0°, 0°); that is, the input axis of the gyro was directed along the m 2 , 
or roll, axis. A comparison, at high gain, between this and the basic system is 
shown in figure 15 . It can be seen that the stability of the system with gyro 
oriented (0°, 0°, 0°) is much poorer than that of the basic system. The high 
frequency, lightly damped mode which occurs primarily in yaw can be anticipated 
from the arguments leading to sketch (e). A similar investigation using low 
gains showed no detectable differences between the behaviors of the two systems - 
again in general agreement with sketch (e) . 


Effect of Control Axes Misalinement 


The effect of a +5° pitch control axis misalinement on the behavior of both 
low and high gain systems is shown in figure l6. It can be seen that with the 
low gain system the pitch inertia-wheel speed and pitch angle reach steady-state 
values in about 50 minutes. On the other hand, with the high gain system there 
is no indication of an approach to the steady-state condition even after 200 min- 
utes. In fact, with this magnitude of pitch axis misalinement, the pitch inertia 
wheel of the high gain system attains its maximum design speed before a steady- 
state condition occurs. Thus, whereas the pitch inertia wheel of the low gain 
system only attains a speed of 80 rpm, that of the high gain system reaches its 
maximum design speed of 1800 rpm in 8 hours and must then be despun back to 
zero - a process which must be repeated at least every 8 hours. 

A possible difficulty may occur with a low gain system because of the angle 
of pitch deviating from zero to a value equal to the misalinement angle. Sup- 
pose, for example, that it is desired to maintain the pitch attitude angle to 
within 3° of zero and that at 3° the coarse control system operates to reduce the 
attitude error. The coarse control system would have to come into operation at 
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least once every 25 minutes to maintain the pitch attitude within the specified 
limits. This could result in a higher rate of energy expenditure by the coarse 
control system than continual inertia-wheel speed desaturation with the high gain 
system. If, however, the specifications are such that the coarse control system 
threshold is greater than the anticipated pitch control axis misalinement , then 
the low gain system may have advantages from the point of view of expenditure of 
energy for control. 

A series of computer runs was made wherein the systems were released from 
various initial conditions. It was observed that, for both low and high gain 
systems, the transients with and without the + 5 ° pitch control axis misalinement 
were almost indistinguishable, thus verifying that small positive pitch control 
axis misalinements have very little influence of stability. 

Shown in figure 17 is the transient behavior when the basic high gain system 
with a pitch control axis misalinement of -5° is released from an initial condi- 
tion of 5 ° of roll. A comparison between figure 17 and figure 15 shows that, 
again, misalinements have very little influence on stability. If, however, the 
gyro orientation is changed from (0°, -6°, 0°) to (0°, -2°, 0°) then (see 
fig. 17) the -5° pitch control axis misalinement causes a high-frequency very 
lightly damped mode which appears chiefly in the yaw degree of freedom. That 
this loss of stability is due to the misalinement and not to the change of gyro 
orientation can be seen from figure 18 which shows the system behavior in the 
absence of misalinements. This demonstrates the importance of pitch axis mis- 
alinements on stability and indicates the sensitivity of the system to gyro 
orientation. A similar investigation carried out using the low gain system 
showed that pitch axis misalinements, of the order of magnitude considered here, 
have very little influence on the behavior of the system. 

The effect of a + 5 ° roll control axis misalinement on the behavior of both 
low and high gain systems is shown in figure 19. It can be seen that a roll axis 
misalinement is not so important as a pitch misalinement since it results in only 
relatively small variations in the orientation angles and inertia- wheel speeds. 
Computer runs comparing the initial transients with and without roll axis mis- 
alinements demonstrated that roll axis misalinements have a negligible influence 
on stability. 


Yaw Channel Failure 

It has been shown previously that, if the yaw channel fails, the system can 
be stabilized by the introduction of a constant speed pitch inertia wheel. It 
was also suggested that, provided the required angular momentum of the constant 
speed wheel is sufficiently small, biasing the speed of the pitch control inertia 
wheel could provide a solution. The feasibility of a scheme of this kind has 
been investigated for the satellite under consideration here. A value of h 2 ^ 
equal to 0.2 ft -lb -sec was selected. This is exactly equal to one half of the 
assumed pitch control wheel capability (appendix D) and also satisfies the 
expression 



h gK ^ 

WoW 



h sK ^ 


= 1 
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thus assuring stability. The root locus diagram for the roll -yaw stability, 
using this value of h^, is shown in figure 20. It can be seen that with the 
low gain system (i.e., Kp = 0.3 ft-lb-sec/radian) the stability is relatively 
good, the damping factor for both modes of oscillation being of the order of 
O.35. It should be noted that this is better than the damping factor of the low 
frequency mode of the fully operating system (damping factor of 0.13). The tran- 
sients following initial conditions of 5° of roll and 0.005°/sec of roll rate are 
shown in figures 21 and 22. It can be seen that a yaw channel failure makes the 
system more sensitive to disturbances which influence roll rate but it operates 
satisfactorily and offers a method of extending the life of the satellite in the 
event of a yaw channel failure. 

The stability of the high gain system (i.e., K<p = 7^.4 ft-lb-sec/radian) is 
poor, as might be expected. The root locus diagram of figure 20 shows that not 
only is the low-frequency mode, for all practical purposes, undamped but there 
is also a negative real root very close to zero (s' = -0.00255) • The implication 
is that any disturbance to either the roll or yaw degrees of freedom will cause 
the yaw angle to oscillate with orbital frequency (and zero damping) while the 
mean angle of yaw tends to drift away from zero. This conclusion is verified by 
figures 23 and 24 which show the transients following initial conditions of 5° 
angle of roll or 0.005°/sec rate of roll. Thus this scheme when applied to a 
high gain system cannot be regarded as providing an adequate solution to the yaw 
channel failure problem. 


Yaw Bias Response 


The response of both the low gain and high gain systems to a yaw bias 
command is shown in figure 25- The command signal used is based on a Nimbus 
satellite requirement represented here as a change of yaw angle from +37 • 5° to 
-37. 5° in about 10 minutes. With the type of time linear command signal shown 
in figure 25, the performance of the high gain system is undoubtedly superior. 

The change of angle of yaw from +37-5° to -37-5° takes about 12 minutes for the 
high gain system as compared with about 37 minutes for the low gain system. 

Because of the coupling between the degrees of freedom, an angle of yaw 
produces steady-state angles of roll and pitch as discussed in earlier sections 
of this report. This effect is particularly evident with the low gain system 
where, as can be seen from figure 25, the angle of roll changes by about 8 due 
to the yaw bias maneuver, that is, from +4° to -4°. The effect on both roll and 
pitch is predicted with good accuracy by equations (71) and (72) . With the high 
gain system the steady-state roll and pitch angles are reduced, relative to those 
of the low gain system, by a factor of 255 (i.e. , ratio of high gain to low gain; 
see eqs . (71) and (72)) and are therefore too small to be detected in figure 25. 

If, in the case of the satellite considered here, roll deviations of the 
order of 4° are tolerable during the yaw bias phase, then the yaw bias require- 
ment can be met by reshaping the command signal. For example, if the command 
signal is changed to that shown in figure 26, then the low gain system meets the 
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requirements. The feasibility of this approach depends, of course, on the 
stringency of the yaw bias requirement which may, in some cases, be sufficiently 
severe to preclude the use of a truly low gain system. 


SUMMARY OF RESULTS 


A dominating feature of the dynamical behavior of the earth-oriented 
satellites considered is that the earth pointing axis rotates with orbital 
angular velocity with respect to inertial space, which results in the roll and 
yaw degrees of freedom being strongly coupled. This coupling, in general, pre- 
cludes an analysis on a unidimensional basis. One important exception is when 
the condition (ig - Ii) - (h 2 j[/co 0 ) = 0 is satisfied. In this case the stability 
in roll and yaw are independent. An investigation of this case reveals many 
fundamental aspects of the dynamical behavior of earth-oriented satellites. 

With the type of control system considered here, there exists a set of gains 
for maximum damping (termed the low gain system) . An increase in the gains above 
those for maximum damping causes a progressive reduction of the damping of a mode 
of oscillation with a frequency close to orbital. This mode is closely related 
to the stability of the vehicle angular momentum, which in turn is dictated pri- 
marily by the gravity torques acting on the vehicle. In fact the magnitude of 
the gravity torques largely determines the maximum rate at which angular momentum 
can be dumped into the orbit. 

The stability of the system is independent of whether or not the gyro input 
axis is located out of the satellite roll-yaw plane. A rate gyro with a small 
time constant oriented, relative to the roll axis, at a small angle above the 
satellite roll pitch plane provides the best method of controlling yaw. An inte- 
grating gyro can be used to control yaw and with the most favorable orientation 
results in a low gain system with a performance comparable with that using a rate 
gyro. Even with the most favorable orientation a high gain system using an inte- 
grating gyro has poorer stability than one with a rate gyro. One problem with 
the integrating gyro is that the most favorable orientation from the point of 
view of stability results in a reduced sensitivity to angle of yaw. 

Control axes misalinements are more serious in pitch than in either roll or 
yaw. A pitch misalinement can cause instability if the rate gyro input axis is 
oriented too close to the roll axis. In addition, pitch misalinements may cause 
excessively high rates of expenditure of stored energy by the coarse control 
system for both low and high gain systems , the former because of large deviations 
in angle of pitch, the latter through the need for repeated inertia-wheel despin- 
ning. The low gain system could show advantages if the anticipated pitch control 
axis misalinement angle were less than the satellite coarse control system pitch 
attitude threshold limit. 

Yaw bias requirements are most easily met by the high gain system, but the 
low gain system in many cases can be made to operate satisfactorily by suitable 
shaping of the yaw command signal. One problem with the low gain system is that 
large steady-state angles of yaw produce steady-state angles of pitch and roll. 
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A constant speed pitch, wheel will stabilize the system when the yaw channel 
fails and is particularly satisfactory when the system gains are low. The system 
is still theoretically stable when high gains are used, but it is so poorly 
damped that it cannot be regarded as satisfactory. A constant speed pitch wheel 
will not stabilize the system when the roll channel fails. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 13, 1963 
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APPENDIX A 


DERIVATION OF THE GYRO TRANSFER FUNCTION 

The expression for the gyro gimbal angle 0g (eq. (2)) may he rewritten in 
the form 

IgRg + Go9g = -IgCCsi^l + C22^2 + ^ 23 ^ 3 ) 

+ H r [(Cii - C3x0g)(Ji + (Cx2 - C320g)(^2 + ( C x3 - C330g)(d3] 

+ AMg - K0 g (Al) 


where K is the torque generator feedback constant (K = 0 for a gyro operating 
in the integrating mode, AMg is any spurious torque acting about the gyro gimbal 
axis and Cii, Cis, etc . , 1 are related to C11, C12, etc . , by the f ollowing 
expressions 


C21 - C 2 1 - t d/m C22 + 0 d/m C 23 

C22 = C 2 2 + ^ d / m C21 " C P d / m C23 

C23 = C 2 3 - %/m^ 21 + ^d/m^ 22 

C'n = On - t d/m Ci2 + 0 d/m Ci3 I 

C 12 = C 12 + t d/m Cn - cp d/m Ci3 

C i 3 = C 13 - 0 d / m Cxx + cpd/m c i 2 


(A2) 


Note that if there are no misalinements, C11 = On, etc . 


In the case of small satellite "body angles and rates, the relationship 
between the inertial angular velocities and Euler angles and rates as expressed 
by equations (3) becomes, to a first approximation. 


“1 = %m/o - “oV/0 

w 2 = ^m/o - w o / 

w 3 = V/o + ^oTm/o 


(A 3 ) 


■••Note that the quantities G X2 y etc., are the direction cosines 

relating the c d axes to the m d axes. 
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Substituting equation (A3) into equation (Al) , rearranging the terms, and 
deleting those involving products of small quantities results in the following 
expression: 

Ig0g + ^D^g + (K - H r o) 2 oC32)% = H r [(C ii^m/o + ^ 13C0 ° C f ) m/ o ^ + ^ 12 ^®m/o " w o) 

+ ( c i3t m/0 - Cii w 0 t m / 0 ) ] + AM g (AU) 

where the terms in Ig on the right side of equation (Al) have been omitted 
because H r is numerically very much greater than Ig (the ratio is usually of 
the order of 10 3 ) . 

Equation (A4), written in terms of the Laplace transform variable s, becomes 

[igS 2 + Cps + (K - Hj>U)qC32 ) ] 0 g = H r [(CixS + Ci3U) 0 )cp m/ . 0 + C' 12 s 0 m / o 

+ (C13S - Ci 1 Wo) 1 l r m / 0 ] + AMg - H r Ci2 w O (A- 5 ) 

Since, in general, Ig/Cp « Cp/(K-H r u> 0 C32) the expression IgS 2 +Cps+(K-H r u) 0 C32) 
can be written, with good accuracy, in the form 

(K - IW^) ♦ l) ( K _ -^c - + 0 


The time constant I g /Cp is usually extremely small, that is, of the order of 
0.001 second. It can therefore be assumed that the term [(l g s/Cp) + 1] does not 
significantly affect the system dynamics. Equation (A5) can then be written in 
the following form: 


B„ = 


Hr 


K -H r W 0 0 32 


AMo 


(CiiS+Ci3Uo)cp m /o+Cx2se m / 0 +(Ci3S-Ci 1 w 0 )t nl/0 + -jj— 'C 12 W 0 


T_S + 1 


(A 6) 


where 

t C D 

Tg _ K - H r to 0 C 3 2 


It can be seen from equation (A 6) that the only difference between the 
dynamic behavior of gyros operating in the rate or the integrating mode arises 
through differences in the time constant T g . As might be expected, the differ- 
ences in Tg can be very large. For example, a typical gyro has a value of T g 
in the integrating mode of about 2000 seconds, while in the rate mode it can be 
made, say, 0.2 second or even smaller, if desired, simply by choosing the appro- 
priate value for the torque generator feedback constant K. The value of Tg in 
the integrating mode can be changed by changing the value of the damping constant 
Cp, although this may involve major design changes in the gyro. 


It should be noted here that the gyro used in the integrating mode is stable 
only if H r (j0 o C 3 2 is negative. Since co 0 is positive and H r is positive, this 
implies that C32 must be negative. This means that the positive direction of 
the gyro spin reference axis must make an angle of less than 90° with respect to 
the negative d 2 direction of the gyro reference frame. This condition is not 
generally necessary in the case of the rate gyro since K is normally much 
greater than H r u) 0 C3 2 . 
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APPENDIX B 


ALTERNATIVE DERIVATION OF THE CHARACTERISTIC EQUATION 


In the section of the report dealing with the stability of the linearized 
equations, the characteristic equation is derived by the usual method of taking 
the determinant of the coefficients. The disadvantage of this approach is that 
it is a purely formalistic procedure which, in some cases, can obscure some of 
the physical facts related to the final form of the characteristic equation. It 
is proposed here to use an alternative method of derivation to reveal some of 
these physi'cal implications . 


The set of linearized homogeneous equations is given below in terms of the 
Laplace transform variable s' appropriate to a time scale r = u> Q t. 


IlS ' 2 - k(l 3 - I 2 ) - 


n sK 


-(ii - i 2 + i 3 )s' . ?ksl ? + 52 


W Q J GdO 


[I 2 s» 2 + 3(ii - i 3 )]e 


s’ ho 


di - I 2 + I 3 )s f + ^ 


-Kq 3 Fcp(d)oS t )d> 


I 3 S' 2 + (I 2 - II) 


^2K 


¥ + 


0 )q 


- K ^ FqCwoS 1 )0 


K^F^(a) 0 s’ )(C 11 s’ + 0^3)$ + K\jfF^(w 0 s' )C 12 © 


+ K\jrF^(a) 0 s’ )(ci3s’ - Ci 1 )Y 


_ 0 
w o 

(Bl) 
= 0 
(B2) 

+ i!lsw = 0 
0 ) o 

(B3) 
= 0 
(B4) 

= 0 

(B5) 

+ = 0 
“° (B 6 ) 


where $>, ©, and ¥ are the Laplace transforms of cp m/0 , ^m/o* f ‘°BB°' w ' s 

immediately from equations (B2)and (B5) that 


I 2 s - 2 + 3(li - Is) + 


K 0 F@(w o s' )s' 


u. 


0=0 


(B7) 


Consider the two quantities £ and t] defined by the following equations 


+ + 

da - II) - 77 s 

0) o 

L w o J 

+ h3w + 

r ( i 3 . l2) + hd 

Wq 

L w o J 


¥ 


(B 8 ) 

(B9) 
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The quantities t, and t] are, in reality, the components of angular momentum 
about the bx and b 3 axes, respectively. Substituting equations (b 8) and (B9) 
into equations (Bl) and (B3) yields 


s'£ - q - 3(ls - I 2 )« = 0 


(BIO) 


s’t) - £ 


= 0 


(Bll) 


Substituting for 0 from equation (BIO) into equations (B4) and (B8) and for 
h lw /oj 0 from equation (B4) into equation (B8) yields 


.■«{ - s-, + m . 3jj 

I iUq 


(la - Ii) - 


h aK 
Wq _ 


¥ = 0 (B12) 


Substituting for £ from equation (Bll) into equation (B12) 

..2 


; '(S' 2 + 1 + 3Ji) + ^ (M ° s ') (s ' 

Ill0 o 


T1 + 3Jl 


da-n) 


U r 


¥ = 0 


(B13) 


Substituting for <E> from equation (BIO) into equations (B9) and (B 6), and for 
h 3w /(jJ 0 from equation (b 6) into equation (B9) 


t\ + 


(Is - la) + 


h 2K 


Wf 


♦ CU) 


o J 3JiI 


3JxI] 


+ K^(u 0 s')CU© - [I 3 s* - K|-F\|r(o) 0 s ')(Cx3 s 1 - Cix)]¥ = 0 (Bl4) 

Substituting for £ from equation (Bll) into equation (Bl4) and rearranging 


- — 13Ji - 
3Jx 


(l 3 - I 2 ) + 


~ >a 1} + K^cOoS’Hc'xis' + C' 13 ) ( — 2 - 1} 




u 0 J I z 


- K^F^OOqS’ )Ci20 + I 3 


K\|rF\j/-(to 0 s' )(Ci 3 s' - Cj_x) 


Is 


Wow if (I 2 - Ii) - (h 2 K/ u, o) = equation (B13) reduces to 


>’ (s' 2 + 1 + 3J X ) + 


K cp F C p(( J o 0 s')(s' 2 + 1)' 
I 1 W 0 


T ] = 0 


1=0 

(B15) 


(B16) 
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thus demonstrating that with [(l 2 - Ii) - (h 2 k/ w o )1 = 0 " the equation 


s'(s' 2 + 1 + 3Ji) 


^(coqsMCs' 2 + 1) 

I i«*>o 


(B17) 


governs the stability of the yaw component of angular momentum r\ and, by equa- 
tions (Bll) and (BIO), the stability of the roll component of angular momentum (; 
and the roll angle $. Thus equation (B17) provides all the necessary informa- 
tion for deciding whether or not the system is capable of "dumping" its angular 
momentum (roll -yaw) into the orbit. 

A disturbance introduced through the angle of yaw Y cannot influence © 
nor, if (I2 - Ii) - (h 2 K/ w 0 ) = can ^ influence t], since both of these varia- 
bles are uniquely determined by equations (B7) and (Bl6). It follows that the 
stability in yaw is governed only by the term in Y in equation (B15) ; that is, 


K^(co 0 s')(Ci 3 s - Cji) 
Is 


(B18) 


This argument also shows that since C' 12 only occurs in equation (B14) it can 
have no influence on the stability of the system. However, while the terms in 
© and T] in equation (B15) do not influence the stability of the system, they do 
have an influence on the transient motion in yaw; that is, a disturbance in 
either 0 or q produces a disturbance in Y. 
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APPENDIX C 


SOME PROPERTIES OF THE DENOMINATOR POLYNOMIAL IN EQUATION (59) 


Consider the polynomial 


s 


t 4 


+ S 


r 2 


1 + 3Ji + 





+ 





= 0 


when 


J x > 0 


To ease the subsequent discussion, the polynomial will he written 


s' 4 +Bs ,2 +C=0 


where 


B 


1 + 3Ji + 





C 



h gK N 

w 0^3/ 




(Cl) 


(C2) 

(C3) 


It is well known that the most stable set of roots 1 which can be possessed 
by a polynomial of the form given in equation (Cl) is when all the roots are 
imaginary. The necessary and sufficient conditions for this are: 


B > 0 
C > 0 
B 2 - 4C > 0 


Equation (C3) can be written in the following form: 


h 2K 


h sK 


U0I3/ WqI 3 


J 3 + 4j + (IU 1 J 3 - C) = 0 


(CU) 
(C5) 
(c 6) 


(C7) 


1 The term "most stable set of roots" is defined to mean the set of roots 

with the least upper bound for the real parts. 
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and has real roots provided 



+ 4J i 



- 4 ^ (4J x J 3 - C) > 0 

13 


or 



4 c > 0 

±3 


(C8) 


Thus if condition (C5) is satisfied, then equation (C7) always has two real roots 
(i.e., 1 and (hgj(;/io 0 l 3 ) 2 ) . These roots clearly satisfy the following 

conditions : 




4J X 




and 




< toi 






(C9) 


(CIO) 


Equation (C2) can he written in the following form 


B = 1 + 3Ji 


1 



h gK ^ 

w 0 l3y. 


+ C 


(Cll) 


Thus if a value of h 2 £ is selected which satisfies condition (C9) and, there- 
fore condition (C5)j» then it follows from equation (Cll) that 

B > 1 + C > 0 (C12) 

so that condition (c4) is satisfied. Furthermore it follows from condi- 
tions (C12) and (C5) that 
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B 2 - 4C > B 2 - 4(B - 1) = (B - 2) 2 > 0 


(C13) 


so that condition (c6) is satisfied. 

The solution of equation (Cl) in terms of s' 2 is 


. t 2 _ 


-B ± n/b 2 - 4C 


Considering first the positive value of the square root and using condition (C12) 


( s '2) = -B + n/b 2 - 4C > -B +n/b 2 - 4(b - 1) 

1 o — o 


= -1 


Thus 


Six = -Sx2 = Kj 


where K < 1 


(C14) 


The remaining roots are derived "by taking the negative value for the square root, 
that is , 


(s , 2) = : B -JEZm < -B_- JE lj = _( B _ 1} < _ c 

\ 2 — 2 — 


Thus 


S 21 = “ s 22 ~ ^ j <s/~C 
where L > 1 


(C15) 


To sum up, therefore, for every value of which satisfies the 

condition 


[hoth (I 2 - Ii) and 4(l 2 - I 3 ) ] 

0) o 


all the roots of equation (Cl) lie on the imaginary axis. Furthermore the roots 
are symmetrically disposed relative to the real axis and the two on the positive 
imaginary axis are located so that one is always equal to or less than j and 
the other always greater than or equal to j >J~C. 
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APPENDIX D 


SATELLITE BASIC DATA 
Inertias 


Principal inertia in roll I]_ 
Principal inertia in pitch. I2 
Principal inertia in yaw I3 


190 slug -ft" 
146 slug -ft" 


120 slug -ft" 


Inertia Parameters 


* ■ M 

I 3 

0.137 

J2= I 2 

I 3 

0.460 

’ l 2 l3 

II 

- 0.367 


Control -System Transfer Functions 


Fcp(S) = 
F 0 (s) = 
*\|r(s) = 


1 + T cp s 


(1 + aT cp s)(l + T m s) 

1 + T e s 

(1 + CGT 0 s)(l + T m s) 

1 

(1 + TgS ) ( 1 + T m s) 
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where the time constants, etc., have the following values: 


Motor time constant T m 

38.4 sec 

Gyro time constant t (rate mode) 

0 

0.2 sec 

Gyro time constant Tg (integrating mode) 

2000 sec 

Compensation network time constant = t q 

7 sec 

Compensation network lead-lag constant a 

0.1 

Additional Motor Flywheel Characteristics 

Maximum flywheel speed 

1800 rpm 

Stall torque 

0.0104 ft 

Maximum momentum 

0.4 ft -lb 


Gyro Orientation Matrix 


Since for a gyro working in the integrating mode, it is essential that C 32 
he negative (see appendix A). The c-^ and d-^ reference frames are initially 
alined so that the spin reference axis of the gyro lies along the negative d^ 
axis, that is, initially, 


Cl" 


0 

rH 


X 

C2 

= 

0 0 1 


d 2 

C3 


0-10 


_d 3 


The Euler sequence orienting the Cj_ frame relative to the <±j_ frame is pitch, 
yaw, roll, that is, 0 C /^, ^c/d' Tc/d' This yields the following relationships 

for those elements of the [C] matrix which enter into the gyro equations: 


Cn = cos i|r c/d cos 0 c/d 

C12 = Sin \|r c / d 

C 13 = sin 0 c/d cos \|r c /a 

C 3 2 = -cos t c/d cos cp c/d 

It is assumed throughout that the orbit is circular and the orbital rate 
u) 0 is 0.001 radian/sec. This is equivalent to a satellite altitude of about 
580 statute miles. 
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TABLE I.- COMPARISON BETWEEN APPROXIMATE AND EXACT ROOTS OF ROLL-YAW CHARACTERISTIC EQUATION 


Low gain Kp = O .3 ft -lb -sec/radian, = 300 ft -lb -sec 2 / radian 


Principal 
degree of 
freedom 

Approximation (l) 
Derived from equa- 
tions ( 36 ) and ( 37 ) 

Approximation (2) 
Equations ( 36 ) and (37) with 
true open loop imaginary poles 

Exact 


Real 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

Roll 

- 0.1009 

- 1.062 

-0.1242 

-1.077 

-0.1376 

- 1.052 

Roll 

-.1009 

1.062 

- . 1242 

1.077 

-.1376 

1.052 

Roll 

-24.61 


-24.61 


-24.76 


Roll 

-1430 


-1430 


-1430 


Roll 

-1.469 


-1.422 


-1.468 

-.9164 

Yaw 

-2.133 


-2.133 


-1.468 

.9164 

Yaw 

-30.76 


-30.76 


- 31.20 


Yaw 

-4993 


-4993 


-4993 



High gain 

Kqj = 76.4 ft-lb-sec/radian. 

= 76,400 ft-lb-sec 2 /radian 


Roll 

-.0005105 

- 1.000 

-.0006210 

- 1.000 

-.0005162 

- 1.000 

Roll 

-.0005105 

1.000 

-.0006210 

1.000 

-.0005162 

1.000 

Roll 

-48.17 

-93.26 

-48.17 

-93.26 

-46.90 

-93-43 

Roll 

-48.17 

93.26 

-48.17 

93-26 

-46.90 

93.43 

Roll 

-1360 


-1360 


-1360 


Yaw 

-2508 

-1588 

-2508 

-1588 

-2509 

-1593 

Yaw 

-2508 

1588 

-2508 

1588 

-2509 

1593 

Yaw 

-9-461 


-9.461 


-9.490 



Satellite orbit 



o 3 is along a line from the vehicle center of mass to the center of 
the earth. 

o 2 is perpendicular to the orbital plane (measured positive in the 
reverse direction to the orbital angular momentum vector). 

o is perpendicular to o 3 and o 2 and completes the right handed 
axes system. 

Figure 1.- Orbital reference frame. 
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Figure 2.- Roll stability. 
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Figure 4 .- Pitch stability 
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Figure 5 .- Response from initial condition 

orientation (O c 


cp = 5°j rate gyro (t~ = 0.2 sec) with 
, -6°, 0°). 
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Low-gain system 
K^=K 0 = 0.3 ft— lb — sec/radian, 
K^=300 ft-lb-sec 2 /radian 


High -gain system 
K<£ = K 0 = 76.4 ft-lb-sec/radian 
K^ = 76,400 ft -lb-sec 2 /radian 
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Figure 6. - 
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Figure 7* “ Response from initial condition 0 = 5°i rate gyro (t„ = 0.2 sec) with 

orientation (0°, -6°, 0°). 
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K£ =K 0 = I .Oft-lb- sec/radian, 
K + = I OOOft-lb-sec 2 /radian 


<^,deg 
0,deg 
'/'.deg 
, deg/sec 
9 , deg/sec 
(//,deg/sec 

w w,- rpm 

w w 2 * rpm 

"w 3 - rpm 



Figure 8.- Response from initial condition cp = 5°; rate gyro (T g = 0.2 sec) with 

orientation (0°, - 6 °, 0°). 
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Integrating — gyro (T g = 2000 sec) 
Orientation (0 o ,-6°,0°) 



Integrating —gyro (r g = 2000 sec) 
Orientation (0°,-70 o ,0°) 


-H h— 12.5 min 



Figure 

Kcp 


9. - Response from initial condition cp = 5 °; low-gain system 
= Kq =0.3 ft -lb -sec/radian, = 300 ft -lb -sec 2 / radian. 
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Figure 10. - Response from initial condition 
Kq> = Kg = 0-3 ft -lb -sec/radian, = 
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- Response from initial condition cp = 5°; high-gain system 
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Figure 12.- Response from initial condition ^ = 0.005°/sec; high -gain system 
Kip = Kg = j6.k ft -lh -sec/radian, = 76,400 ft -lb -sec 2 / radian. 
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Figure 14.- Response from initial condition t = 0.02°/sec; high-gain system 
IQp = Kq = 76.4 ft-lb-sec/radian, Ity = 76,400 ft -lb -sec 2 / radian. 
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Figure 15*- Response from initial condition cp = 5°; high -gain system 
= Kg = 76.4 ft -lb -sec/radian, Ify = 76,400 ft-lb-sec 2 /radian. 
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Figure 16. - Effect of 5° axis misalinement in pitch (initial orientation angles 
and rates zero); rate gyro (t c = 0.2 sec) with orientation (0°, -6°, 0°). 
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Figure 17- - Response from initial condition 9 = 5°; high-gain system 
K<p = Kg = 76.4 ft -lb -sec /radian, = 76,400 ft-lb-sec 2 /radian. 
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Figure l8.- Response from initial condition cp = 5°; high -gain system 
Kip = Kg = 76.4 ft -lb -sec /radian, = 76,400 ft-lb-sec 2 /radian. 
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- Effect of 5° axis misalinement in roll (initial orientation angles and 
zero); rate gyro (Tg = 0.2 sec) with orientation ( 0 °, “ 6 °, 0 °). 
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Figure 20.- Roll -yaw stability after yaw channel failure (hgj ^ = 0.2 ft-lb-sec) 
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21.- Response from initial condition cp = 5°; low-gain system 
= 0.3 ft-lh-sec/radian, = 300 ft -lb -sec 2 / radian; 
gyro (Tg = 0.2 sec) with orientation (0°, -6°, 0°). 
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Figure 23-- Response from initial condition cp = 5°; high -gain system 
K^p = Kg = "j6.b ft-lb-sec/radian, Ify = 76,400 ft-lb-sec 2 /radian; 
rate gyro (Tg = 0.2 sec) with orientation (0°, -6°, 0°). 
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Figure 24.- Response from initial condition 9 = 0.005°/sec; high-gain system 
Kq;, = K@ = 76.4 ft-lh-sec/radian, = 7^,400 ft-lh-sec 2 /radian; rate gyro 
(Tg = 0.2 sec) with orientation ( 0 °, - 6 °, 0 °). 
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Figure 25-- Response to yaw lias signal (rate 
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Figure 26.- Response to yaw bias signal (low-gain system 

Kq, = Kg = 0-3 ft-lb -sec/radian, = 300 ft-lb-sec 2 /radian) ; 
rate-gyro orientation (0°, -6°, 0°). 
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